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ABSTRACT 

We consider some of the spectral and polarimetric signatures of the gravitational confined detona- 
tion scenario for Type la supernova explosions. In this model, material produced by an off-center 
deflagration (which itself fails to produce the explosion) forms a metal-rich atmosphere above the 
white dwarf surface. Using hydrodynamical simulations, we show that this atmosphere is compressed 
and accelerated during the subsequent interaction with the supernova ejecta. This leads ultimately 
to the formation of a high- velocity pancake of metal-rich material that is geometrically detached from 
the bulk of the ejecta. When observed at the epochs near maximum light, this absorbing pancake pro- 
duces a highly blueshifted and polarized calcium IR triplet absorption feature similar to that observed 
in several Type la supernovae. We discuss the orientation effects present in our model and contrast 
them to those expected in other supernova explosion models. We propose that a large sample of 
spectropolarimetric observations can be used to critically evaluate the different theoretical scenarios. 

Subject headings: hydrodynamics - supernovae: general - polarization 



1. INTRODUCTION 

After decades of study, the mechanism of the explosion 
of a white dwarf in a Type la supern ovae remains uncer- 
tain (Hill ebrandt fc Niemeverll20Q0lL Several of the pro- 
posed theoretical mo dels (|ArnettHl969l: iKhokhlovl 2001; 
iReinecke et al.ll2002l) have failed to produce objects with 
the energetics and chemical comp osition compatible with 
observations ( Hof lich et al.l2002) . These two characteris- 
tics are best captured by the so-called delayed- detonation 
mode ls (JKhokhlov 199l HGamezo et all2004HPlewa et alJ 
2004). In such models, a mild ignition occurs near the 
center of the accreting, massive white dwarf, and sparks 
a deflagration (subsonic flame). During the subsequent 
evolution, a deflagration to detonation (supersonic reac- 
tive wave) transition (DDT) takes place. 

Despite the promise of the DDT models, it remains 

unclear how a transition to detonation occurs - in all 

standard DDT models calculated so far, the detona- 

^ ^^^^^^^ 

tion w as triggered artificially. However, IPlewa et alJ 
( 2004) have proposed a gravitational confined detonation 
(GCD) model in which a detonation naturally follows a 
slightly off-center ignition. In this scenario, the deflagra- 
tion takes the form of a single bubble, buoyantly rising 
to the stellar surface. At the bubble's breakout, the sur- 
face layers of the star are laterally accelerated and begin 
sweeping across the stellar surface, converging opposite 
the breakout point. The subsequent compression of the 
colliding streams and thermalization of the kinetic energy 
triggers a detonation. The authors speculate that, as in 
the standard DDT case, the detonation will consume the 
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remaining fuel, producing an energetic explosion with 
chemically stratified ejecta. The stability of the GCD is 
currently the subject of more detailed numerical study. 

One of the most striking properties of the GCD model 
is that the products of the deflagration are brought to 
the surface of the white dwarf prior to detonation. This 
material constitutes a small fraction of the stellar mass, 
and is rich in metals. This compositional "pollution" of 
the outer stellar layers is reminiscent of a peculiar fea- 
ture noticed in the spectra of some Type la SNe. In a 
handful of objects, observers have identif ied a highly- 
blueshifted Ca II IR triplet absorption (lHatano et al.| 
1999; |Li et al.| iWang et al.| EH |Gerardy et al.| 

2004), indicating absorbing material moving at veloci- 
ties rsj 20,000 km s _1 , much higher than that charac- 
teristic of other spectral lines. Given that velocity is 
proportional to radius in expanding supernova atmo- 
spheres, the high-velocity (HV) calcium absorption in- 
dicates a component of absorbing material geometrically 
detached from the bulk of the ejecta. Spectropolarimetry 
of SN 200 lei further showed that the HV absorption was 
highly polarized, indicating that the absorbing material 
was distri buted aspherically with perhap s a clump-like 
geometry (Ei et alJBflflj |Kasen et alj|Il|). 

In this letter we demonstrate that the ejecta structure 
characteristic of the GCD scenario can naturally explain 
observations of the peculiar HV calcium absorption fea- 
ture in Type la supernovae. 

2. METHODS 

We studied the post-detonation evolution of the GCD 
model using a simplified numerical setup in which spher- 
ical supernova ejecta runs into an aspherical metal-rich 
atmosphere. Fo r the supernova eject a structure, we used 
the W7 model (|Nomoto et al.lll984l) with an initial ra- 
dius of 3 x 10 8 cm. Around the ejecta, we placed an 
ellipsoidal metal-rich extended atmosphere representing 
the bubble of burned material expelled during the GCD 
breakout. The composition of this material was taken 
to be the oxygen-burned composition #4 from Table 3 
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Fig. 1. — Hydrodynamical simulation of the model supernova ejecta interacting with the extended atmosphere. Panels (a)-(c) show 
the density in log scale; panel (d) depicts the magnitude of the velocity. The density scale is shown by the colorbar in panel (c) while 
the velocity scale is shown by the colorbar in panel (d). The black contour line corresponds to a calcium number abundance of 0.01. (a) 
Density distribution at the beginning of the simulation (t = s). (b) Density distribution at t = 0.3 s. Notice the remarkable deformation 
of the calcium-rich material in the upper part of the computational domain, (c) Density distribution at the final time (t = 1.24 s). Notice 
that the calcium-rich region has been strongly compressed into a pancake-like structure, (d): Velocity magnitude at the final time. The 
calcium-rich absorber is seen moving at velocity ~ 21, 000 km s — 1 with a substantial velocity gradient across the structure. 



of Khok hlov et alJ (|1993fi which consists of 57% silicon, 
27% sulfur, 7.1% iron, 2.7% calcium and small amounts 
of other elements. The extended atmosphere had an axis 
ratio of 1.2, semi- major axis length of 6 x 10 8 cm and 
was centered at 3 x 10 8 cm. The density of the at- 
mosphere decreased exponentially with a scale length of 
3 x 10 8 cm. The atmosphere was uniformly expanding 
about its center with velocity proportional to radius and 
reaching 5000 km s _1 at the outer edge. The ejecta and 
ellipsoidal extended atmosphere were embedded in an 
ambient medium of pure helium of density 10 -4 g cm -3 . 
The initial model was isothermal with T — 10 T K. This 
overall configuration closely repr e sents that seen in the 
GCD calculations of iPlewa et all (|2004fl . 

The hydrodynamic evolution was calcul ated with the 
adap tive mesh refinement code flash ((Frvxell et alJ 
2000). The initial model was defined on a 2-D cylin- 
drical grid covering the region up to 2.56 x 10 10 cm in 
radius and from —2.56 x 10 10 cm to 2.56 x 10 10 cm in the 
z-direction. The ejecta was centered at the origin of the 
grid. The maximum resolution of the simulation, used 
only to resolve strong flow structures in the dense regions 
of the model, was equal to 62.5 km. We used a reflecting 
boundary condition at the symmetry axis and allowed 
for free outflow otherwise. We use d a Helmholtz equa - 
tion of state, an isol3 composition (jFrvxell et alJ l2000). 
and a multipole solver with 10 terms in the expansion to 
account for self-gravity. 

Spectra of the hydrodynamical models were calculated 
using a multi-dimensional Monte Carlo radiative trans- 
fer code (|ThomaJ2003 lKasenll20"ol . The opacities used 
in the calculation were electron scattering and bound- 
bound line transitions. Ionization and excitation were 
computed assuming local thermodynamic equilibrium, 
where the temperature structure of the atmosphere was 
determined self-consistent ly usi ng an iterat ive approach 
enforcing radiative-equilibrium (|Lucvl 119*9?^ . Line inter- 
actions were treated in the Sobolev approximation and 
included both absorption and scattering according to a 
equivalent two level atom scheme with thermalization pa- 



rameter e = 0.01. Monte Carlo photon packets were ini- 
tially emitted from a spherical inner boundary surface 
located at v = 6000 km s _1 and according to a black- 
body distribution with T = 9000 K. The packets were 
initially unpolarized but acquired polarization by elec- 
tron scattering. Line scattered light was assumed to be 
unpolarized due to complete redistribution. To construct 
the emergent spectra and polarization, escaping photon 
packets were collected into one of 100 angular bins. 

3. RESULTS 

We have calculated several models assuming different 
mass of the extended atmosphere. Here we discuss the 
evolution for one case in detail (m atm — 0.008 M ) and 
describe the differences in the other cases as appropriate. 

The initial setup of our model is shown in Fig. GJa). 
The W7 ejecta occupies the innermost region of the grid, 
shown in shades of red. The yellow ellipsoidal region sur- 
rounding the ejecta represents material expelled during 
the deflagration. The black contour lines enclose the re- 
gions where the calcium number abundance exceeds 1%. 
In particular, the contour in the top portion of the image 
marks calcium produced in the deflagration; the inner 
calcium-rich ring is that of the W7 nucleosynthesis. 

The overall hydrodynamic evolution in our model is 
best described as ejecta impacting a stratified atmo- 
sphere, with shock-related effects playing a minor role. 
In Fig. nfb) (t = 0.3 s), a mostly spherical expansion 
of the supernova shock can be seen in the lower part 
of the computational domain. In the upper part of the 
domain, the ejecta impacts and deforms the extended 
atmosphere. By t = 1.24 s (FigH^c)), that impact has 
lead to the compression of the calcium-rich material into 
a pancake-like structure of radius ~ 30,000 km, located 
near z ~ 2 x 10 9 cm. By that time, the ambient medium 
surrounding the atmosphere has been completely over- 
run by the forward supernova shock. The reverse shock 
can be seen as the highly contrasted spherical structure 
starting at (r, z) = (0, —2.6 x 10 9 ) cm. 

Beginning at t « 1, the overall expansion becomes in- 
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Fig. 2. — Comparison of the synthetic model spectra to Type la 
supernova observations. Top panel: Spectral observati ons of two 
Type la SNe near maximum light; SN 2001el (red line. lWang et all 
2003) shows a strong HV calcium abs orption at 8000 A while 
SN 1994D (black line, Patat et al. 1996) does not. Bottom panel: 
Synthetic model spectra at 20 days. A HV calcium absorption is 
clearly seen when looking straight down on the calcium-rich pan- 
cake (viewing angle = 0°, red line), whereas none is seen from 
the opposite side (0 = 180°, black line). 
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Fig. 3. — Synthetic polarization spectra of the model at 20 days. 
By convention, a positive (negative) polarization level signifies po- 
larization aligned parallel (perpendicular) to the symmetry axis. 
The black lines show the polarization level while the gray lines 
show the corresponding (arbitrarily scaled) flux spectrum. The in- 
sets illustrate how the HV calcium line polarization is caused by 
the partial obscuration of the supernova photosphere (light blue 
disk) by the calcium-rich pancake (dark blue object). 



creasingly homologous, and we stop our calculation at 
t = 1.24 s. Fig. IHd) shows the total velocity in the 
upper portion of the computational domain at the final 
time. The calcium-rich pancake structure shows a sig- 
nificant velocity gradient across its body (Av/c « 0.02). 
The central region of the pancake (r = 0) moves with 
velocity spanning 17,000-24,000 km s _1 , while near the 
outer rim the velocity is higher (22,000-28,000 km s _1 ). 

In the case of a more massive atmospheres, the over- 
all evolution and resulting morphology are similar to the 
case described above. The major difference is the final 
velocity of the calcium-rich pancake. For a model with 
^atm = 0.016 Mq, the inner part of the pancake moves 
with velocities 14,000-20,000 kins" 1 , while the outer 
edge moves at 21,000-26,000 km s _1 . For an even more 
massive atmosphere (m atm = 0.08 M ) the velocities 
are still lower, increasing from 10,000-12,000 km s _1 to 
17,000-21,000 km s _1 across the pancake. In this case 
the calcium-rich pancake almost overlays the region of 
intermediate mass elements in the W7 eject a. 

The synthetic spectrum of the model at a time when 
the supernova is near maximum light has been obtained 
after homologously expanding the eject a to 20 days 
(lower panel of Fig. 0). Given the enhanced abundance 
of intermediate mass elements and relatively low temper- 
ature (T « 5500 K), the pancake is opaque in the Ca II 
IR triplet lines. For a viewing angle, in which the ob- 
server looks directly down upon the pancake (0 = 0°), the 
pancake obscures the supernova photosphere, creating 
a broad and highly blueshifted HV calcium absorption 
feature near 8000 A. The HV feature seen in the model 
comp ares well to that observed in SN 200 lei (jWaner et alJ 
2003), shown in the top panel of Fig. [2] For larger view- 
ing angles, the pancake obscures less (or none) of the 
photosphere, and the HV feature is weaker or absent in 



the model spectrum. This orientation effect may explain 
why in some super novae, such as SN 1994D ([Patat et alJ 
1996), a HV calcium feature at maximum light is seen 
only weakly or not at all. 

The aspherical geometry of the calcium-rich pancake 
also leads to significant polarization over the HV calcium 
feature (Fig. EJ). Light in the supernova ejecta becomes 
polarized by electron scattering. By convention, a pos- 
itive (negative) polarization level signifies polarization 
oriented parallel (perpendicular) to the symmetry axis. 
Because the electron-scattering photosphere in our model 
is essentially spherical, the observed polarization cancels 
in the continuum. However, the HV pancake may par- 
tially obscure the underlying photosphere, leading to a 
non-zero polarization over the HV calcium feature. This 
effect is illustrated by the insets to Fig. El 

The line polarization created in this way depends sen- 
sitively upon the viewing angle. For 6 = 0°, the polar- 
ization cancels completely due to the azimuthal symme- 
try (Fig. Ef a ))- As is increased, the absorber selec- 
tively reveals horizontally polarized light from the lower 
rim of the photosphere, resulting in negative line polar- 
ization of order - 1% (Fig. CJb)). For = 45°, the 
model closely reproduc es the HV line po larization peak 
observed in SN 2001el ((Wang et alJl200fl . For larger 0, 
vertically polarized light from the side edges of the pho- 
tosphere dominates, and the line polarization changes 
sign (Fig. 0Jc)). For > 110°, the pancake no longer ob- 
scures the photosphere, and leaves no obvious signature 
in either the flux or polarization spectrum (Fig.Efd)). 

4. DISCUSSION 

We have demonstrated that the observations of a pe- 
culiar HV calcium feature in Type la supernovae are nat- 
urally explained within the GCD model. In the model, 
burned material expelled during the breakout of the de- 
flagrating bubble forms an extended atmosphere above 
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the stellar surface. That atmosphere is subsequently 
compressed and accelerated during the hydrodynamic in- 
teraction with the supernova eject a. This leads to the 
formation of a high-velocity calcium-rich pancake with 
cross-section comparable to that of the underlying pho- 
tosphere. The partial obscuration of the photosphere by 
the pancake results in a HV calcium absorption feature 
with significant line polarization. 

In our model, the blueshift of the absorption fea- 
ture depends sensitively upon the mass of the expelled 
material. For the fiducial case studied here (m a t m = 
0.008 M ), the pancake material spans the velocity 
range 17,000-24,000 km s _1 , and is geometrically de- 
tached from the bulk of the supernova ejecta. This com- 
pares well with the velocities inferred from the HV cal- 
cium feature of SN 2001el. As the atmosphere mass is 
increased, the absorbing pancake moves at lower veloc- 
ity and eventually blends with the region of intermediate 
mass elements in the supernova ejecta. In such a case 
we might expect a very different observable signature, in 
which the pancake material increases the strength and 
blueshift of several of the normal Type la spectral fea- 
tures. This could provide an orientation-dependent ex- 
planation for the unusually high velocities and peculiar 
velocity evolution of the normal featur es in some Type la 
SNe (e.g., SN 1984 A (iBranchl IT987h and SN 2002bo 
(|Benet,t,T i et, a1.ll200lV 

The large size of the absorbing pancake results in the 
distinctive orientation effects present in our model. The 
strength of the HV flux absorption depends only upon 
the fraction of the photosphere obscured by the pancake, 
and hence decreases monotonically with 0. The line po- 
larization, in contrast, is maximal when either a large or 
a small part of the photosphere is covered. Such behavior 
could be tested with a large sample of spectropolarimet- 
ric observations of SNe la; in particular, one could study 
the correlation between the HV calcium polarization level 
with the depth of the flux absorption. 

Although the HV calcium feature is the most profound 
signature of the absorbing pancake, the opacity due to 
numerous iron and titanium lines creates a modest flux 
depression in the wavelength region 3500-4500 A. This 
leads to a variation of the B-band magnitude with view- 
ing angle of order one tenth of a magnitude. As dis- 
cussed in|Branch et al] ^M) and|Thomas e{jp §004), 
the temporal evolution of this opacity should also affect 
the shape of the supernova light curve. Therefore, the 
presence of the calcium-rich absorber is another intrin- 
sic source of SNe la photometric diversity, of possible 
relevance to their cosmological application. 



The presence of a HV calcium absorber can possibly 
be explained in other explosion scenarios. For example, 
the large bubbles of burned material present in the pure 
deflagration and standard DDT models may in princi- 
ple produce a number of HV absorbers. This case can 
be distinguished from the case of single absorber in the 
GCD model using a rich sample of spectropolarimetric 
observations. In particular, our model predicts a much 
smaller fraction of supernovae showing a persistent HV 
calcium feature, as well as the aforementioned polariza- 
tion correlations. Note that multi-epoch observations 
may be necessary to separate the signatures of the HV 
ab sorber from the "transient ionization" effect discussed 
bv lGerardv et all (|2004l) . whereby the recombination of 
Ca III to Ca II in the cold outer layers of ejecta may also 
give rise to a (short-lived) HV IR triplet absorption in 
the epochs prior to maximum light. 

For several SNe la, polarization in the continuum has 
also been detected, indicating a global asymmetry in the 
bulk of the supernova ejecta. In SN 2001el, the con- 
tinuum polarization angle differed from that of the HV 
calcium feature, suggesting that the orientation of the 
HV absorber deviated from that of the ejecta. One pos- 
sible explanation for this difference is that the ejecta ac- 
quired a separate, large-scale asy mmetry due to the in- 
teraction with a companion star ((Marietta et alJ 2000: 
Kas en et al.ll200J) . Another interesting possibility is to 
consider the GCD framework with rotation of the pro- 
genitor included. In such a case, the trajectory of the 
buoyantly rising deflagrating bubble is not expected to 
be aligned with the rotation axis of the progenitor. This 
lack of correlation follows from the fact that the con- 
vective core of the white dwarf is expected to produce 
ignition seed points in a largely stochastic manner. Such 
a scenario should be studied in the future by means of in- 
tegrated multi-dimensional hydrodynamical simulations. 
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